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ABSTRACT 

We present a detailed study of the X-ray energy and power spectral properties 
of the neutron star transient IGR J17191-2821. We discovered four instances of pairs 
of simultaneous kilohertz quasi-periodic oscillations (kHz QPOs). The frequency dif- 
ference between these kHz QPOs is between 315 Hz and 362 Hz. We also report on 
the detection of five thermonuclear type-I X-ray bursts and the discovery of burst 
oscillations at ^ 294 Hz during three of them. Finally, we report on a faint and short 
outburst precursor, which occurred about two months before the main outburst. Our 
results on the broadband spectral and variability properties allow us to firmly establish 
the atoll source nature of IGR J17191-2821. 

Key words: Keywords: accretion, accretion disks — binaries: close — stars: individ- 
ual (IGR J17191-2821, XTE J1747-274) — stars: neutron — X-rays: stars 



1 INTRODUCTION 

Neutron star low-mass X-ray binaries (NS-LMXBs) have 
been extensively observed with the Rossi X-ray Timing Ex- 
plorer (RXTE) during the last 13 years. These observations 
have led to important discoveries, such as persistent and 
intermittent pulsations in accretion-powered millisecond X- 
ray pulsars, nearly coherent oscillations during X-ray bursts 
and strong quasi-periodic variability on millisecond time- 
scales (the so called kilohertz quasi-periodic oscillations; kHz 
QPOs). 

The kHz QPOs are relatively narrow peaks in the power 
spectrum. They sometimes occur in pairs and are thought to 
reflect motion of matter around the neutron star at the inner 
edge of the accretion disk (see, e.g. iMiller et al.lll998l ). Up 
to date, kHz QPOs with similar characteristics have been 
detected in about 30 neutron star sources. While it is often 
assumed that at least one of the QPOs reflects the Keplerian 
mot ion at the inner edge of the accretion disk (see e.g. review 
by Ivan der "k lis 2006), models have to still satisfactorily 
explain the presence and characteristics of both QPOs. 

Direct detection of coherent or nearly coherent pulsa- 
tions is the only available method to measure the neutron 
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star spin period in LMXBs. Te n NS-LMXBs out of more 
than 100 known (|Liu et al.ll2007l ) have shown coherent mil- 
lisecond pulsations in their persistent emission; these sys- 
tems are known as accretion-powered millisecond X-ray pul- 
sars (AMXPs; see Wii nands . 20051 . for a r eview of the first six 
AMXPs discovered; for the last four s ee lKaaret et al. 20061. 
Krimm et al ]|2007l . lCaseH"a" et al ]|2008l and lAltamirano et al.l 



2008al ). Sixteen sources to date (including the one pre- 
sented in this paper) have shown nearly coherent millisec- 
ond oscillations during t hermonuclear Type-I X-ray bursts 
I Watts et al.ll2008l . I2OO9I '). As the X-ray burst evolves, the 
oscillation frequency typically increases by a few Hz, ap- 
proaching an asymptotic value (ubo) which is stable for a 
given source from burst to burst. This asymptotic frequency 
is thought to trace within a few Hz the spin frequency {vs) 
of the neutron sta r (IStrohmaver et aLlll99d '). The AMXPs 
SAX J1808-3658 j Chakrab artv et'aLri2003h. XTE J1814 



338 (|Strohmaver et a l. 2003), Aql X-1 jCasella et al 



2008 ) 



and most recently HETE J 1900. 1-2455 (| Watts et all 120091 ) 
have all shown that Vs — vbo, strongly supporting the idea 
that Vbo is a good tracer of Vs. 

Since the discovery of burst oscillations and two simul- 
taneous kHz QPOs in neutron star LMXBs, it was sug- 
gested that there was a relation between the neutron star 
spin frequency, Vs, and the kHz QPO frequency difference 
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Figure 1. Swift/XRT image of IGR J17191-2821. The black 
circle represents the INTEGRAL error for the source position. 
The source showed an short (days) outburst precursor or flare on 
March 8. The long outburst studied in this work occurred in May 
(Swift/XRT detection is showed as bright source on the top-left 
of this figure). 



Av = Vu — vi- Although the frequency of the QPOs was 
found to vary (150 <iyi< 900 Hz and 350 <yu< 1200 Hz), 
initial measurements revealed that Av was consistent with 
being constanlQ and e qual to the asymptotic burst oscil- 
lation vbo (see, e.g., Istrohmaver et all Il996l ). This was 
the main motivati on for beat-frequen cy models such as the 
sonic-point model (|Miller et al.lll99"8l ). which proposed that 
Vu reflects the Keplerian frequency at the inner edge of the 
disk, and that ut was the beat between and Vs- As new 
observations revealed more sources showing both vbo and 
twin kHz QPOs, the r elati on became more complex (see 
iMendez fc BeUonill2007l and Ivan der Klisll2008l for a detailed 
historical overview). As a summary, NS-LMXBs are often 
classified as fast or slow rotators, depending on whether the 
spin fr equency is higher or lower than ~ 400 Hz, respec- 
tively (|Muno et al.ll200lh . It was found that generally the 
fast rotators follow Av ~ '^s/2i while slow rotators f ollow 
Av ~ Vs (see IWiinands et al.ll2003l : lLinares et al.ll2005l , and 
references therein). This of course implies that there is no 
one-to-one relation between Av and Vg. Furthermore, it is 
now known that Av is not always consistent with being con- 
stant for sources with known vbo (see, e.g ., iMendez et al.l 
1 19981 : IJonker et aLll2002l : iBarret et allbood . and references 
within) . 

The question of whether Av and Vs are physically 
related is still under debate. While current data might 
still be compatible with a bimodal relation, recent results 
suggest that this migh t not be the case. For example, 
IMendez fc Bellonil l|2007 ^ suggested that Av and Vg are un- 
related and that the division between fast and slow rotators 
may be just an effect of the lo w number of so urces show- 
ing both phenomena (s ee also lYin et al] |2007| ) . Recently, 
Istrohmaver et al] (|2008l ) reported the discovery of burst os- 
cillations at 414.7 Hz in the LMXB 4U 0614+091; ff con- 



^ At least in those sources in which vbo was known. For example, 
Ai/ in the neutron star L MXB Sco X-1 is known to vary (see, e.g. 
Ivan der Klis et al.|[l997l) . 



firmed, these results would imply a spin frequency that is 
inconsistent with either Av ~ Vs/2 or Av ~ Vs relations. 



1.1 IGR J17191-2821 

IGR J17191-2821 was discovered by the IBIS/ISGRI instru- 
ment aboard INTEG RAL during observatio ns of the galac- 
tic bulge monitoring (iKuulkers et al.ll2007bl) c onducted be- 
tween March 2 and 4. 2007 jTurler et al.ll2007l ). The source 
was detected in the 20-40 and the 40-80 keV bands with sig- 
nificances of 10.4 and 4(t, and fluxes of 8.9±0.9 and 5.8± 1.4 
mCrab, respectively. The position of the source was first re- 
ported as (RA,DEC) = (259°. 77, -28°. 35) (J2000) with an 
accuracy of 2.5 arcmin. 

On March 3 , 2007, RXTE Galactic bulge scans 



(|Swank fc Markwardtll200ll ) dete cted a source at a p osition 
consistent with that reported by iTurler et al.l (l2007h at an 
inten sity of 10± 1 mCrab in the 2-10 keV band (jSwank et al.l 
I2OO7I ). Follow-up RXTE observations performed three to 
four days later did not detect the source with a 3(t upper 
limit of 1.2 mCrab. Eight years of Galactic bulge scans were 
reanalyzed for contributions from a source at the position 
of IGR J17191-2821, but no flares brighter than 2 mCrab 
were found. These non-detecti ons showed that th is source is 
active relatively infrequently l| Swank et al.ll2007l ). 

During further Galactic bulge scan observations on 
April 29 and May 2, 2007, IGR J17191-2821 was detected 
again, but at a level of 30 and 70 mCrab (2-10 keV). 
This suggested that the previous detection was a fiare or 
outburst precursor. Several Astronomical Telegrams at this 
time communicated on the discovery of Type I X -ray bursts 
dKlein-Woh et al."2007b': 'Markwa rdt et al.ll2007l ). burst os- 
cillations (Markwardt ct al. 2007.) and several episodes of 
kHz QPOs (|Klein-Woh et allbOOT^ ). 

In this work, we present an intensive study of Type I X- 
ray burst characteristics, burst oscillations and kHz QPOs of 
this newly discovered neutron star LMXB IGR J17191-2821. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Light curves, color diagrams and timing 
analysis 

We use data from the Rossi X-ray Timing Explorer (RXTE) 
Proporti onal Counter Array (PGA; for instrume nt informa- 
tion see IZhang et all Il993l : Ijahoda et al] I2OO6I ). To study 
the long-term (days/months) Lx behavior of the source, 
we use d the PGA monitoring ob servations of the galactic 
bulge l|Swank fc Markwardtl [20011 ) . These observations are 
performed nine months of the year (as parts of the months 
of November, December, January and June are excluded due 
to solar constraints). The accuracy in the position of the 
PGA bulge scans is about 15 arcmin; the light curves are 
given in the ~ 2 — 10 keV energy band. 

To study the short-term (minutes or less) variability, we 
use PGA pointed observations. For IGR J17191-2821 there 
were 18 observations in one data set (92052-10) containing 
~ 2.5 to ~ 10 ksec of useful data per observation. We use 
the 16-s time-resolution Standard 2 mode data to calculate 
X-ray colors. Hard and soft color are defined as the 9.7-16.0 
keV / 6.0-9.7 keV and 3.5-6.0 keV / 2.0-3.5 keV count rate 
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ratio, respectively, and intensity as the 2.0-16.0 keV count 
rate. We removed Type I X-ray bursts from the data as 
well as corrected by deadtime effects and for the contribu- 
tion of the background. We normaliz ed colors and intensi- 
ties by those of the Crab Nebula fsee iKuulkers et aL 1994 
van Straaten et all l2003l . see table 2 in lAltamirano et al.l 
2008bl for average colors of the Crab Nebula per PCU). 

For the Fourier timing analysis we used the Event 
mode E_125us_64M_0_ls. Leahy-normalized power spectra 
were constructed using data segments of 128 seconds and 
1/8192 s time bins such that the lowest available frequency 
is 1/128 Hz and the Nyquist frequency 4096 Hz. No back- 
ground or deadtime corrections were performed prior to the 
calculation of the power spectra. Unless stated explicitly, in 
our fits we only include those Lorentzians for which we can 
measure the integrated power with an accuracy of at least 
3(T, based on the (negative) error in the power integrated 
from to CX3. For the kHz QPOs, we report the centroid 
frequency vo, the full width at half maximum (FWHM) and 
the rms amplitude. The quoted errors use Ax^ = 1.0. The 
upper limits quoted in this paper correspond to a 95% con- 
fidence level (Ax^ = 2.7). 



trum, which means that for a purely random Poisson pro- 
cess, the powers follow a x^ distribution with 2 degrees of 
freedom (IStrohmaver fc M arkwardt"l999^. We searched the 
30-4000 Hz frequency range in the 2-60 keV band and in 
narrower energy bands (Sec. I3.2|l . 

We also created energy spectra every 0.25 sec from the 
Event mode (E_125us_64M_0_ls) data of all the PCUs that 
were on during the burst. Given the high count rates dur- 
ing the peak of the bursts, we corrected each energy spec- 
trum for dead-time using the methods suggested by the 
RXTE teanJ3. For each energy spectrum, we created the 
corresponding response matrix using the latest information 
available on the response of the instrument at the relevant 
times. As is common practice, we used as background the 
energy spectrum of the persistent emission taken seconds be- 
fore each burst. (We used 100 sec of the persistent emission 
to calculate the spectrum. However, we found no significant 
differences in the fits when the persistent-emission before or 
after the burst was used, or when using data-segments of 
different lengths -between 100 and 500 seconds-). We used 
a black-body model to fit the resulting burst spectr£0. 



2.2 Energy spectra of the persistent emission 

For the PCA, we used the Standard 2 data of PCU 2, which 
was active in all observations. The background was esti- 
mated using PCABACKEST version 6.0 (see FTOOLS). 
We calculated the PCU 2 response matrix for each obser- 
vation using the FTOOLS routine PCARSP VlO.l. For the 
HEXTE instrument, spectra were accumulated for cluster B 
(as cluster A stopped rocking in October 2006), excluding 
the damaged detector and averaging both rocking directions 
to measure the background spectrum. Dead time corrections 
of both source and background spectra were performed us- 
ing HXTDEAD V6.0. The response matrices were created 
using HXTRSP V3.1. Both for PCA and HEXTE, we fil- 
tered out data recorded during, and up to 30 minutes after 
passage through the South Atlantic Anomaly (SAA). We 
only used data when the pointing offset from the source was 
less than 0.02 degrees and the elevation of the source respect 
to the E arth was grea ter than 10 degrees. Using XSPEC 
V11.3.2i (lArnaudlll9"9^ '). we fitted simuhaneously the PCA 
and HEXTE energy spectra using the 3.0-25.0 keV and 
20.0-200.0 keV energy bands, respectively. We used a model 
consisting of a disk blackbody and a power law, absorbed 
with a n equivalent Hydrogen column density of 0.3x10^^ 
cm (Hi em- Woh et al.ir2007ah . which gave a good fit in all 
the observations {y^ /dof < 1.1). 



3 RESULTS 

3.1 Position of the source 

Swift observed the source twice on March 8, 2007 (at 02:37 
UT for a total of ~ 800 sec, and at 10:44 UT for a total 
of ~ 1600 sec - both observations where performed in PC 
mode). These are the only Swift observations performed be- 
fore the bright outburst (see Section [3. 4p . The observations 
were taken in Photon Counting mode and did not show a 
source within the INTEGRAL error circle with an upper 
limit of 0.0021 cuts s~^ (at a 95% confidence level). As- 
suming a galactic absorption Nh of 3.4 x lO^'^ cm~^ and a 
photon index of 1.8, this countrate translates into an upper 
limit of ~ 8.7 x 10~^* erg s"'^ cm~^ on the unabsorbed flux 
(approximately 0.004 mCrab in the 2-10 keV range). We 
found a source in the first Swift observation (i.e. on March 
8 at 2:37 UT) located at (RA, DEC) = (259.8114, -28.3005) 
(J2000, with an error of ~ 9 arcsec), i.e., at a distance of 

^ http: / /heasarc. gsfc.nasa.gov/docs/xte/recipes/pca_deadtime. html 
^ We assumed that the X-ray spectra after the persistent emis- 
sion has been subtracted are Planckian and that the observed 
luminosity of the source is: 

L = A-KoT^B?, 

so the unabsorbed bolometric X-ray flux may be determined using 



2.3 Type I X-ray bursts 

We examined the Standard 1 mode data (2-60 keV, 0.125 
seconds time resolution, no energy resolution) of the 18 ob- 
servations for Type I X-ray bursts; we found 5 episodes (see 
Table [TJ . We searche d each burst for coh erent pulsations 
using the statistic l|Strohmaver fc Markwardt. 19991 1 . We 
computed Z\ (i.e. assuming that the signal is sinusoidal) 
throughout the bursts using a sliding 2 seconds window 
with a step of 0.125 sec. The statistic has the same 
statistical properties as a Leahy normalized power spec- 



F^,l=aT^{R/Df, 

where a is the Stefan-Boltzmann constant, T is the black-body 
temperature, R the neutron star photosphere radius, and D the 
distance to the source. The ratio [R/D)^ is the normalization 
of the black-body model we used (bbodyrad - see Xspec manual 
for details). We note that X-ray burst spectra are generally well 
described by black-body emission, however, the emission from the 
neutron star and its environment (e.g. accretion disk) is expected 
to be m ore complex than simple black-body emission (see, e.g., 

Ivan Pa radiis 1982; L ondon et al 1984 : Kuulkcrs et al...200a . and 

references therein). 
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Table 1. X-ray bursts in IGR J17191-2821 



Number 


ObsID 


Start time of the 


Num. of 


Flux" 


Osc. 


Osc. rms amplitude'' 






burst (UT, 2007) 


PCUs on 


(10~* erg cm~^) 


(yes/no) 


(2-17 keV) 


1 


92052-10-01-00 


May 4 02:32:06 


2 


2.5 ±0.3 


yes 


6.9 ± 0.6% 


2 


92052-10-05-00 


May 7 02:39:39 


2 


2.9 ±0.3 


yes 


5.0 ± 1.0% 


3 


92052-10-03-01 


May 7 05:51:12 


2 


2.7 ±0.3 


no 


< 3% 


4 


92052-10-06-00 


May 8 17:08:53 


3 


2.6 ±0.2 


yes 


10.2 ±1.5% 


5 


92052-10-06-01 


May 8 20:38:44 


2 


2.1 ± 0.3 


no 


< 3% 



":Bolomctric peak flux. 

*": Integrated amplitude of the oscillations in the 2-17 keV range. 

about 3.5 arcminutes from the position of the source as mea- 
sured by INTEGRAL. We measured an average count rate of 
0.0087 cuts s~^, which corresponds to an unabsorbed flux of 
3.6 X 10"^^ erg s"^ cm"^ or 0.02 mCrab (2-10 keV, assuming 
the same Nh and photon index as above). In the second Swift 
observation this source is not significantly detected (a total 
of three photons within a 10 pixel source error circle and 
zero from a background region 3 times as large). We place a 
95% confidence upper limit on the 2-10 keV unabsorbed flux 
of ~ 2.0 X 10"^^ erg s"^ cm"^ (employing PIMM^ for an 
absorbed power-law spectrum with A^^ = 3.4 x 10^^ cm~^ 
and a photon index of 1.8, and applying the proscription for 
low number statistics given bv lGehrels, 198 6). In Figure[T]we 
show a Swift/XRT image during IGR J17191-2821 outburst 

(May 10th, 2007)^ 

As noted bv iKlein-Woh et all (|2007al l. the formal IN- 
TEGRAL error circle on the position of IGR J17191-2821 
would suggest that the faint source we detected with Swift 
is unrelated to IGR J17191-2821 (see Figure [T|. However, 
given the systematic uncertainties in the INTEGRAL posi- 
tion, both sources are probably one and the same. We note 
that this is not the flrst case in which the true position of a 
transient laid outsid e the reported INTEG RAL-IBIS /ISGRI 
error circle (see e.g. iKuulkers et al]|2007al ). 



3.2 Thermonuclear X-ray bursts and the distance 
to the source 

We found 5 Type I X-ray bursts (see Table [T]). All bursts 
showed similar temperature, radius and bolometric flux pro- 
files (not shown) . The temperature and fiux profiles were all 
single peaked, and reached their maxima within a second. 
The maximum temperature (kT) was always between 2 and 
3 keV and the peak bolometric fluxes were in the 2 — 3 x 10~* 
erg s~^ cm~'^ range (Table [T}. In all cases, the black-body 
radius remained approximately constant after the peak and 
was usually constrained between 5 and 10 km (assuming a 
distance of 11 kpc). None of the bursts showed indications 
of photospheric radius expansion. By using the highest mea- 
sured bolometric peak flux of 3 x 10~* erg s~^ cm~^ we can 
estimate an upper limit on the distance. We flnd a distance 
D < 11 kpc when using the empirically determined Edding- 
ton luminosity of 3.79 x 10^® e rg s~^ (for bursts show ing pho- 
tospheric radius expansion - iKuulkers et al.ir2003l ). Using a 



^ http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 




Frequency (Hz) 



Figure 2. In the main panel we show the average kHz QPOs of 
all data reported in Table [J] after the shift-and-add methods was 
applied; the lower kHz QPO was shifted to the arbitrary frequency 
of 700 Hz (see Section|33)l . Inset: kHz QPOs in observation 92052- 
10-05-00; no shift-and-add was applied. 

more complex approximatioijfl and standard values for the 
mass and the radius of the neutron star (i.e. Mjvs = 1.4Mq 
and ii = 10 km), we found D < 8.6 and D < 6.6 for hy- 
drogen mass fractions of X = and X = 0.7, respectively. 
Larger values of R give higher upper limits. 

^ The approximation was recently used bv lGallowav et al] (|2008| ) 
to compare a sample of more than a thousand X-ray burst from 
different sources. The distance is given by: 

n — c fi I Fl-^^Bol ~i-l/2 I _MmS_\1/2 

where Mf^g is the mass of the neutron star in solar masses, X is 
the mass fraction of hydrogen in the neutron star atmosphere and 
z(R) is the term that takes into account the gravitational redshift 
at the photosphere (were l + z{R) = (1 -2GMns /Rc'^)~^^'^ , G is 
the gravitational constant, c the s peed of light and R th e radius 
measured at the photosphere — see lGallowav et al.ll200^ . 
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3.3 kHz QPOs 

We searched the averaged power spectrum of each observa- 
tion for the presence of significant kHz QPOs at frequencies 
> 200 Hz. In each case, we fitted the power spectra between 
200 and 4000 Hz with a model consisting of one or two 
Lorentzians and a constant to account for the presence of 
QPOs and Poisson noise, respectively. We found that 12 out 
of the 17 observations show significant QPOs in the 605-1185 
frequency range. In 4 observations we detected 2 simultane- 
ous kHz QPOs (see Tablegl). The lower kHz QPO frequency 
was between 680 and 870 Hz, with single trial significance 
between 5.7 and lOcr. The upper kHz QPO frequency was 
between 1037 and 1085 Hz, with single trial significances be- 
tween 3.0 and 3.7(t. A// showed no significant changes and 
was always consistent with 350 Hz (see Table [2]). 

We tried to better constrai n A^' by using the shift- 
and-add method as described bv lMendez et alT (jl998h . We 
first tried to trace the detected kilohertz QP O using a dy- 
nami cal power spectrum (e.g. see figure 2 in iBerger et al.l 
1 19961 ) to visualize the time evolution of the QPO frequency, 
but the signal was too weak to be detected on timescales 
shorter than the averaged observation. Therefore, for each 
observation we used the fitted averaged frequency to shift 
each lower kHz QPO to the arbitrary frequency of 700 Hz. 
Next, the shifted, aligned, power spectra were averaged. The 
average power spectrum was finally fitted in the range 300- 
2048 Hz so as to exclude the edges, which are distorted due 
to the shifting method. To fit the averaged power spectrum, 
we used a function consisting of two Lorentzians and a con- 
stant to fit the QPO and the Poisson noise, respectively. In 
this case, the averaged Ai/ is 332 ± 16 Hz. In Figure [2] we 
show the shift-and-added power spectrum and a representa- 
tive example of the single observation power spectrum with 
two kHz QPOs (see inset). 



3.4 Outburst evolution 

Figure [3] shows the PCA Ught curve of IGR J17191- 
2821 as seen by the PCA bulge scan monitor program 
llSwank fc Markwardtll200ll ). While a precursor of the out- 
burst (see Section ll.l|l occurred at MJD 54162.6, the 
full X-ray outburst did not start until 54 days later (i.e. 
MJD 54216); it lasted for about 30 days. 

On MJD 54247 (May 27th, 2007) the source was not 
detected anymore with RXTE, and a C handra/HRC-I ob- 
servat ion was performed. As reported bv lChakrabartv et al.l 
iiooj), the source was not detected in the 1.1 ks observa- 
tion within the 30 arcsec of the Swift/XRT position (see 
Section [211}. These authors estimated an upper limit on the 
0.3-10 keV unabsorbed fiux of < 8.3 x 10"^" erg s"^ cm"^. 

In Figure [4] we show the persistent unabsorbed 2- 
200 keV flux (panel a), power law index (panel fe), strength 
of the broad band noise (panel c) and kHz QPO frequency 
(panel d) as a function of time during the ~3 weeks of the 
outburst from which RXTE pointed observations are avail- 
able. The source reached a maximum fiux of ~2.5xl0~^ erg 
cm^ s~^. Assuming D < 11 kpc (see Section [3. 2p . we place 
an upper limit on the outburst peak luminosity of 4x10^^ 
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Figure 3. PCA light curve of ICR J17191-2821 as seen by th e 
PCA bulge scan monitor program llSwank fc Markwardj|2001^ . 
The times of the Swift and Chandra observations are marked with 
vertical lines. Contiguous non-zero measurements are connected 
with lines. 



erg s ^. In panel a of Figure [4] we also plot the (Type I) X- 
ray burst bolometric peak fiuxes at the time they occurred 
(as detected by RXTE). 

Due to the relatively low count rates collected by the 
PCA, the average power spectrum of each observation had 
low statistical quality. However, in the brightest (and soft- 
est) observations we found traces of the so-called very low 
frequency noise (VLFN). As a steep power law rising to- 
wards low frequencies, this VLFN is a typical signature of 
the so call ed "banana branc h" (soft state) of atoll sources 
(see, e.g., Ivan der Klij|20()6l . for a review). When compar- 
ing the results showed in the different panels of Figure |4j we 
found that the 5-50 Hz averaged fractional rms amplitude is 
anti-correlated with the source luminosity, whereas the fre- 
quencies of both kHz QPOs showed no obvious trend. The 
spectral index was clearly anti-correlated with the strength 
of the variability (i.e, correlates with luminosity). This is 
simil ar to what has bee n found in other atoll sources (see, 
e.g., Ivan der Klislliooel . for a review), where the strength 
of the variability and the spectral index trace the changes 
in the timing and spectral state during the outburst. From 
Figure |4j we found that the source was initially in the soft 
(banana) state and showed a failed transition to the hard 
(extreme island) state around MJD 54226. After this, it re- 
brightened and returned to the soft state. On MJD 54228 
(i.e. two days later) the luminosity reached a secondary peak 
and started to decline, while the source gradually transi- 
tioned from the soft (banana) state to the hard (extreme 
island) state. Finally IGR J17191-2821 faded below the de- 
tection limit of RXTE-PCA. The timing and spectral prop- 
erties (and evolution) allow us to firmly establish the atoll 
source nature of IGR J17191-2821. 



3.5 Burst oscillations 

By applying the method (see Section I2.3|l on the 2- 
60 keV X-ray burst data, we discovered highly significant 
nearly-coherent oscillations in 2 of the 5 X-ray bursts (Burst 
1 and 4 in Table [1]). For the remaining 3 bursts, we re- 
peated the Zl analysis using only data in different energy 
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Table 2. Observations and kHz QPOs 





IViJJJ 


Aver. 






Lower 






Upper 








(days) 


Cts/s*^ 


on^ 














Au 






/PCU2 




u (Hz) 


FWHM (Hz) 


rms (%,) 


u (Hz) 


FWHM (Hz) 


rms (%) 


(Hz) 


-01-00 


54224.09 


182 


2 


- 


- 


- 


- 


- 


- 


- 


-02-00 


54225.07 


170 


3-4 


870 ± 1 


11 ± 1.7 


8.5 ±0.4 


1185 ± 50 


220 ± 94 


9.2 ± 1.5 


315 ±50 


-02-01 


54225.14 


188 


3 


866 ± 3 


38 ±6 


9.2 ±0.6 


- 


- 


- 


- 


-03-00 


54226.21 


172 


1-3 


881 ± 10 


132 ± 23 


15 ± 1 


- 


- 


- 


- 


-05-00 


54227.10 


142 


2-3 


681 ± 5 


55.2 ±15 


10.3 ±0.9 


1043 ± 10 




8.4 ± 1.1 


362 ± 11 


-03-01 


54227.23 


148 


2 


730 ± 1 


13.2 ± 2.7 


9.7 ± 0.6 


1075 ± 12 




8.5 ± 1.3 


345 ± 12 


-04-00 


54228.01 


156 


3 


891 ± 3 


37.7 ± 7.5 


9.4 ±0.6 










-04-01 


54228.20 


151 


2-3 


838 ± 1 


25 ±4 


9.3 ±0.4 










-06-00 


54228.70 


129 


2-3 


946 ± 20 




10.6 ± 1.2 










-06-01 


54228.90 


130 


2-3 


702 ± 3 


37 ±7 


10.4 ±0.6 


1037 ± 15 


»S_27 


8.8 ± 1.1 


335 ± 15 


-07-00 


54229.06 


105 


1-2 


884 ± 13 


108131 


15 ±2 










-07-01 


54229.98 


72 


2 


706 ± 10 


83 ±25 


18 ±2 










-05-01 


54230.04 


70 


2 


682 ± 20 




22.4 ± 2.8 










-08-00 


54231.09 


46 


2 
















-08-01 


54231.75 


27 


2-3 
















-09-00 


54232.01 


23 


3 
















-09-01 


54232.08 


21 


3-4 

















Background and deadtime corrected averaged count rate for PCU2; this PCU was on during all observations. 
In case the number of active PCUs changed, we report the minimum and maximum number. 



sub-bands. We found that the oscillations were significantly 
detected also in burst number 2, but only in the 10-25 keV 
range. The oscillations are not significantly detected in ei- 
ther burst 3 or 5. Fractional rms amplitudes (averaged over 
the period the signal was significantly detected) and upper 
limits are given in Table [T] In Figure [S] we show the dy- 
namical power spectra of bursts 1 and 4 (upper and lower 
panels, respectively) in the 2-60 keV range. As can be seen, 
the frequency of the oscillations drifts from ~291Hz (burst 
1) and ~ 292.5Hz (burst 2), to a frequency between 294 and 
294.5 Hz. Maximum Zi powers of ~ 63 were found in both 
cases. This type of frequency drifts are typical for burst os- 
cillat ions during Type I X-ray bursts (see, e.g., IStrohmaveij 
I2OOII . for a review) . 

3.6 The energy dependence of the burst 
oscillations 

We analyzed the energy dependence of bursts 1 and 4 which 
showed the strongest oscillations - see Table [T] We split the 
data in six energy bands (from 2.5 keV to 17 keV), and 
calculated the amplitude of the pulsations in each band. To 
create a pulse profile in each energy band, we selected all the 
data with a significant pulse detection in the corresponding 
power spectrum, and explored the P — P — P space around a 
given initial guess value for the pulse period (obtained from 
the power spectrum) . P and P were initially set to zero. We 
note that P and P do not represent true spin changes of the 
neutron star, but comprise all the frequency variations due, 
primarily, to the burst oscillation drift seen in the data. P 
and P are therefore useful to align the phases of the pulsa- 
tions (folded in a profile of 32 bins) for each energy band. We 
then fitted the pulse profiles with two sinusoids representing 
the fundamental and the 1st overtone of the burst oscilla- 
tions. The errors on the fractional amplitudes are calculated 
using a Ax^ = 1- Upper limits are at a 95% confidence level 



(i.e. using Ax^ = 2.7). We found no significant second har- 
monic in any of the energy bands we chose for either burst, 
with rms amplitude upper limits in the 2.5 — 17 keV range of 
3.7% and 2.0%, for bursts 1 and 4 respectively. The energy 
dependence of the burst oscillation rms amplitude is shown 
in Figure [S] The fractional amplitude of the fundamental 
is consistent with being constant. No significant phase lags 
were detected. 



4 DISCUSSION 

We present an intensive study of the X-ray variability of 
the newly discovered neutron star LMXB IGR J17191-2821. 
Our results allow us to firmly establish the atoll source na- 
ture of IGR J17191-2821. We detect several episodes of ther- 
monuclear X-ray bursts, some of them showing burst oscil- 
lations that imply a neutron star with spin frequency of 
~ 294 Hz. We also detect several instances of kHz QPOs; 
when two were detected simultaneously, the difference in 
frequency is consistent with being constant (See Table [2]). 
The energy and broadband power spectra of IGR J17191- 
2821 evolved in a manner consistent with that seen in other 
neutron star LMXBs: it is soft when the fiux is high and 
hard when the flux is low. Near the end of the outburst, 
IGR J17191-2821 shows the hardest spectra and strongest 
variability, with an rms amplitude above 15%. 

Two months before the main outburst, the source exhib- 
ited a very brief (only days) event which was nearly an order 
or magnitude less luminous. It is unclear what the relation 
of this event is with respect to the main outburst. However, 
we note that similar precursors have been seen before (see, 
e.g.. lDegenaar fc Wiinand3l2009D . To our knowledge, no sys- 
tematic search has been performed in order to quantify how 
common these precursors are and how they can be explained 



IGR J 1719 1-2821 7 




§.1000 
o 
900 

o 

O 800 



54222 54224 54226 54228 54230 54232 
Time (MJD) 

Figure 4. Spectral and timing properties along the outburst of 
IGR J17191-2821. From top to bottom, a): 2-200 keV unabsorbed 
persistent flux (flUed circles) and bolometric peak flux of the type 
I X-ray bursts detected by the PCA (open triangles; peak burst 
flux values are divided by ten for plotting purposes), b): power law 
index obtained from the broad band spectral fits; c): fractional 
rms amplitude of the variability between 5 and 50 Hz in the ~2.5- 
45 keV energy band and d): frequencies of the upper (squares) and 
lower (circles) kHz QPOs detected during the outburst. 



in the commonly used disk-inst ability mode ls proposed for 
outbursts of X-ray binaries (see lLasotall200ll . for a review). 

There are currently 12 LMXBs with reported Av mea- 
surements whose spin can be estimated from either pulsa- 
tions i n their persistent X -ray emission or from burst oscilla- 
tions (Ivan der Klisll2008l . for a recent overview). In Figure [7] 
(up per p anel) we plot Av /vs vs. fs for all these sources 
(cf. Ivan d er Klis 2006; Mcn dez fc Bellonill2007l : fvan der Klis 
bOOa ). The dashed line represents a step function: Av jvs ~ 1 
for the slow rotators (< 400 Hz) and Avjvs ~ 0.5 for 
the fast rotators (> 400 Hz). Although most of the data 
seem to be consistent with this scheme (and our data of 
IGR J17191-2821 are as well), certainly there are points that 
do not follow this relation. A clear example is given by the 
neutron star 4U 0614-1-09. A tentative ~ 414.7 Hz burst 
oscillation frequency wa s recently reported for this source 
IStrohmaver et al]|2008l , note that we quote this value as 





Figure 5. Lightcurves and dynamical power spectra (2-60 keV) 
for bursts 1 and 4 (upper and lower panel, respectively — see 
also Table [TJ . The dynamical power spectra were created using 
overlapping 2 s windows, with new windows starting at 0.125 s 
intervals. We used a Nyquist frequency of 4096 Hz. The contours 
show Leahy normalized powers between 20 and 60, increasing in 
steps of 5. 



tentative since it has been detected only once, the signal 
showed no frequency drift as expected from burst oscilla- 
tions, and it was the first and only detection to date of 
burst oscillations with the Burst Alert Telescope on board 
the Swift telescope). As shown in Figure [7]this burst oscilla- 
tion frequency is very close to the discontinuity of the step 
function (although given the present data, the discontinuity 
of the step function could be anywhere between ~ 360 and 
401 Hz; furthermore, to our knowledge none of the models 
predict the exact value at which Au/nUs should switch from 
1 to 0.5). At the same ti me, Au/vs seems to cover almost the 
complete 0.5-1. range (|van Straaten et"alll2000l ). although 
we note that 93% of the Av/us measurements are in the 
0.69-0.85 range and the average Av/va using all mea sure- 
me nts is 0.767 ± 0.006 (see al so, iBoutelier etaLll2009l ). 

iMendez fc Bellonil |20o3) have recently suggested that 
Au and Va are unrelated and that the division between 
fast and slow rotators may be just an effect of the low 
number of sou rces showing both phenomena (see also 
lYin et al] |2007| ) . In the lower panel of Figure [7] we show 
Av and Vs for the same data plotted in the upper panel. 
The dashed line in this panel corresponds to the aver- 
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Figure 6. Fractional rms amplitude versus energy for the oscil- 
lations in X-ray bursts 1 and 4 (See Table[TJ. Arrows mark upper 
limits at 95% confidence level. 



age Af = 308 Hz iMendez fc Bellonil (|2007| ). As can be 
seen, the range of most sources overlaps with this con- 
stant value, except for the two AMXPs SAX J1808.4-3658 
and XTE J 1807-294, for which Av falls clearly below 300 
Hz. The discrepancy is solved, if the data for these two 
sources ar e multiplied by a factor o f 1.5. This was first sug- 
gested by Mendez fc Bellonil (12007) based on the w orks of 
Ivan Straaten et al.l (120051 ) and 'Linar es et all (|2005l ) Q. We 
note that not a ll AMXPs ar c affec ted b y the same multi- 
plicative factor (|van Straaten et al.ll2005l ). and furthermore, 
that the factor might be independent of whet her the neutron 
star pulsates or not l|Altamirano et al1l2005h . 

Clearly, the present data are not enough to draw any 
final conclusion on the relation between Av and Vs', the de- 
tection of both the spin frequency and Af in other sources 
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Figure 7. In the upper panel we plot the ratio between the in- 
dividual measurements of Av = v-u, — and the spin frequency 
Vs as a function of Vs (c f. Ivan der Klidliooel : iMendez fc Bellonil 
l2007l: Ivan der Klij lioosh . Spin frequencies are estimated from 
the persistent pulsations observed in AMXPs (marked with the 
arrows) or from burst oscillations in the non-pulsating sources. 
The lower panel shows the same data as above, but we plot i\v 
vs. Vs. The dashed lines show the step function S(i's) = 1 for 
Vs ^ 400 Hz; S(i's) = 0.5 for Vs > 400 Hz. The continuous line s 
correspond to a constant Au = 308 Hz llMendez fc Bellonill2007^ . 
For 4U 0614 +09 we used the tentati ve burst oscillation frequency 
of 414.7 Hz l lStrohmay er et al .1200^ and the kHz QPOs measure- 



van Straaten et al.l ||2003). For Aql X-1 we used 
Casella et al .1 120081. and r eferences with i n) an d 



ments reported by 

Us = 550 Hz (e.g. 

the tentative Au measurements reported bv lBarret et al. I l l2008t) 
IGR J17191-2821 data are from this work 



° Previous works have shown that the frequencies of the vari- 
ability components observed in atoll sources follow a univer- 
sal scheme of corre lations when plotted versus Uu (see e.g. 
van Straaten et al ]2P02, 200a: lReig et al.l2004l : rAltamirano et al.l 
2008bl . and references therein). The two AMXPs SAX J1808.4- 
3658 and XTE J1807-294 show similar relations. However, the 
relations for these two sources are shifted with respect to those o f 
the other sources llvan Straaten et al]|2005l : iLinares et al.ll2005h . 
This shift is between the frequencies of the low-frequency com- 
ponents and Uu by a factor around 1.5 (1.45 and 1.59 for 
SAX J1808.4-3658 and XTE J1807-294, respectively) and be- 
tween U£ and Uu by a similar factor. 



is necessary. The new instrument ASTROSAT (an Indian 
multiwavelength Astronomy Satellite), which is planned to 
be launched in 2010, will play a major role in solving this 
issue as it is likely to increase the sample of sources with 
both spin frequency and kHz QPOs measurements. 
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